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We investigate the interaction between three rubidium atoms in highly excited (58p) Rydberg 
states lying along a common axis and calculate the potential energy surfaces (PES) between the 
three atoms. We find that three-body long-range potential wells exist in some of these surfaces, 
indicating the existence of very extended bound states that we label macrotrimers. We calculate 
the lowest vibrational eigenmodes and the resulting energy levels and show that the corresponding 
vibrational periods are rapid enough to be detected spectroscopically. 



PACS numbers: 03.65. Sq, 31.50.Df, 32.80. Ee, 34.20.Cf 
I. INTRODUCTION 



Ultracold Rydberg systems are a particularly inter- 
esting avenue of study. Translationally, the atoms are 
very slow, yet their internal energies are very high. The 
large excitation of a single electron leads to exaggerated 
atomic properties, such as long lifetimes, large cross sec- 
tions, and very large polarizabilities [JJ, which can lead to 
strong interactions between Rydberg atoms jH [3] . Such 
interactions have led to various applications in quantum 
information processes over the past decade (see [4] for a 
comprehensive review) . 

Another active area of research with Rydberg atoms 
is in the area of long-range "exotic molecules" . Such 
examples include the trilobite and butterfly states, so- 
called because of the resemblence of their respective wave 
functions to these creatures. First predicted in [5J, these 
quantum states were detected more recently in [Bj. Also 
of interest are the formation and detection of macroscopic 
Rydberg molecules. In [7], it was first predicted that 
weakly bound macrodimers could be formed from the in- 
duced Van der Waals interactions of two Rydberg atoms. 
However, we have shown more recently [8l [9] that larger, 
more stable dimers can be formed via the strong ^-mixing 
of various Rydberg states. Recent measurements [TO] 
have shown signatures of such macrodimers using an ul- 
tracold sample of cesium Rydberg atoms. 

More recently, the focus of study has moved toward 
few-body interactions, such as between atom-diatom in- 
teractions [TTHT3] and diatom-diatom interactions [TU 
ITS] , Coinciding with this shift, there have been pro- 
posals [16419) for many-body long-range interactions in- 
volving Rydberg atoms. However, these works focus on 
the interactions between one Rydberg atom and ground 
state atoms or molecules. In this paper, we describe the 
long-range interactions between three Rydberg atoms ar- 
ranged along a common axis and provide calculations, 
which predict the existence of bound trimcr states. Here, 
we also present the lowest vibrational energies of these 
bound states, calculated via the oscillation eigenmodes 
of the bound system. 



II. THREE BODY INTERACTIONS 

In [5] and [5] , we predicted the existence of long-range 
rubidium Rydberg dimers by analyzing potential energy 
curves corresponding to the interaction energies between 
the two Rydberg atoms. In these works, we diagonal- 
ized an interaction Hamiltonian consisting of the long- 
range Rydberg-Rydberg interaction energy and atomic 
fine structure in the Hund's case (c) basis set. Each 
molecular state in the basis was symmetrized with re- 
spect to the -Doo/i symmetry of the homonuclear dimer. 

In general, adding a third atom to the interaction pic- 
ture will change the physical symmetry of the system. 
However, to simplify our calculations, we assume that 
three identical Rydberg atoms lie along a common (z-) 
axis (see Fig. [TJa)). This preserves the -Doo/i symmetry 
and permits the use of much of the two-body physics 
on the three-body system. We analyzed this symmetry 
state for three Rb 58s atoms and three Rb 58p atoms 
and found that the 58p case exhibited examples of three- 
dimensional wells, indicating that the three atoms are 
bound in a linear chain. 



A. Basis States 

Obtaining properly symmetrized basis functions for 
the three-atom case is very similar to that of the two- 
atom system (see |20j). but much more technically de- 
manding. We construct the molecular wave functions 
from three free Rydberg atoms in respective states 
\ai) = \n 1 J 1 ,j 1 ,m : j 1 ), \a 2 ) = \n 2 ,l2,j2,mj 2 }, and \a 3 ) = 
|a 3 , ^3, j 3 , m^g), where is the principal quantum num- 
ber of atom i, li is the orbital angular momentum of 
atom i, and rrij i is the projection of the total angular 
momentum j i = t\ + Si of atom i onto a quantization 
axis (chosen in the z-direction). As in the two- atom 
case, we assume that the three Rydberg atoms interact 
via long-range dipole-dipole and quadrupole-quadrupole 
couplings. Here, long-range indicates that the distance 
between each Rydberg atom is greater than the Le-Roy 
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Hamiltonian consists of a three-body long-range Ryd- 
berg interaction energy and atomic fine structure, i.e. 
H iut = Vi3_bodv + Hf 8 . Using the wave functions defined 

R 2 by Equation we write the matrix elements of the 

Hamiltonian as the sums of multiple interactions. Each 
>>' matrix element is defined as: 

(ai;a 2 ;a s ,\V 3 ^ hod y\bi;b 2 ;b 3 ) = 

\a 2 ) \a 3 ) - \ a i a j a k Oj, b k , ) (d) 

— ► i,j,k 

(b) i'j'.fc' 

2 X (Oc + PuQa) (&C + Pb®A) , 



where each summation index is over the total number of 
atoms, i.e. from 1 to 3, P is as before, we have defined 



FIG. 1. (Color online) (a) Three Rydberg atoms lie along 
a common z-axis. The distance between atoms 1 and 2 is 
represented by Ri and the distance between atoms 2 and 3 
is represented by i?2- Each atom is in state \a,i), defined in 
the text. Each bound macrotrimer has two eigenmodes of 
oscillation: (b) uj\ and (c) ml. . (See text) 



&a = 



-1 for cyclic permutations 
for anti-cyclic permutations 

and 

for cyclic permutations 
-1 for anti-cyclic permutations 



radius [21] : 

R LR = 2 [(nr^Mi) 172 + (^ 2 |r 2 M 2 ) 1/2 ] , (1) 

such that there is no overlapping of the electron clouds. 
The properly symmetrized long-range three-atom wave- 
functions take the form: 

hi; a 2 ; a 3 ) = 

-yg [(ki)i|a2>2|a3>3 + |a2)i|a 3 )2|ai)3 

+ kshk^Ms) (2) 

-P (\ai)i\a 3 ) 2 \a 2 )3 + |a 2 >i|ai)2|a3>3 

+ |a3)i|a 2 )2|ai>3)] , 

where P — p(— l)* 1+ ' a +^ 3 , w ith p = +1(— 1) for gerade 
(ungerade) molecular states. 

The basis set consists of the Rydberg molecular level 
being probed (e.g. 58p + 58p + 58p) , as well as all nearby 
asymptotes with significant coupling to this level and to 
each other. All of these states are properly symmetrized 
via Equation ^ according to their molecular symmetry 
fl = rrij 1 + rrij 2 + rrij 3 . In this paper, we consider the 
SI = 1/2 symmetry. 



and we have defined \a\ 'a- a k ) = |ai)i|aj)2|dfc)3, etc. 
In the case that \ai;a 2 ;a 3 ) — |&i; ^2; (i-e. along the 
diagonal of the matrix), the matrix element is given by: 

(ax; a 2 ; a 3 \H- mt \ai; a 2 ; a 3 ) = 

(ai; a 2 ; a 3 \V 3 - ho d y \ai; a 2 ; a 3 ) + Ei 23 , (4) 

with E 123 =E X +E 2 +E 3 , where E k = -\(n k - <5 £ J~ 2 
are the atomic Rydberg energies with respective quantum 
defects Se k . 

The long-range assumption assures that these are 
three free atoms interacting via long-range two- 
body potentials. That is, the transition element 
(a^afa^ \V 3 . hody \b^bfb k ^) given in Equation (§ is 
defined as a sum of two-body interactions: 

<a«af >a«|%-bo*|#W> = 

(af\f\V L {R 12 )\b^bf) 
+ (afaf\V L (R^)\bfbf) 

+ (a| 1) 4 3) |^(i?i3)l4 1) ^ ) ) , (5) 



B. Interaction Hamiltonian 

As in the two-atom case, we construct the interac- 
tion picture for the three-atom system by diagonaliz- 
ing an interaction Hamiltonian in the properly sym- 
metrized basis described in the previous subsection. The 



where R a/ 3 is the nuclear separation between atoms a 
and /?, and L — 1(2) for dipolar (quadrupolar) interac- 
tions. Since we are assuming that the three atoms lie 
along a common axis, each two-body interaction term 
(a^ a^\VL(R a p)\b^b^) defines the long-range transi- 
tion element between the two respective Rydberg atoms. 
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Each transition element is given by [91 [5U] : 
<12|Vb(i2)|34) = 
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where j to t = ji +32 +J3+34, m 
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e i = 2£ l + l, ^ = 2j t + 1, and ft£ 



radial matrix element. 



/c 13 /c 24 
R 2L+1 



(6) 



"J 2 



is the 



III. POTENTIAL ENERGY SURFACES 
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FIG. 2. (Color online) Potential energy surface (PES) cor- 
related to the |56p|, |; 58p|, — 4; 60p|, V) asymptotic state. 
This surface is analogous to a repsulsive potential curve for 
the two-body case: As the distance of either the first or last 
atom in the linear chain is increased, the two local atoms re- 
main repulsed. The color scheme denotes the energy values 
given in GHz, with the scale presented to the right of the plot. 



A. General Cases 

We diagonalize the three-body Hamiltonian at succes- 
sive values of R\ and i?2 , resulting in a series of potential 
energy surfaces (PES), where each surface corresponds 
to a different molecular asymptote in the basis. In each 
of the plots shown, R\ represents the distance between 
atom 1 and atom 2, R 2 represents the distance between 
atom 2 and atom 3, both in ao (see Fig. [l|a)) and the 
energy is measured in GHz. The color scheme for the 
energy values is given in the scales to the right of each 
plot. 

As a result of the large ^-mixing that occurs be- 
tween the Rydberg atoms, these surfaces have in- 
teresting topographies. For example, Figures [2] 
and [3] illustrate potential surfaces analogous to two- 
dimensional repulsive and attractive curves, respec- 
tively. The repulsive PES shown in Fig. [2] corre- 
sponds to the |56p|, |; 58p|, — |; 60p|, j) state, while 
the attractive PES shown in Fig. [3] corresponds to the 
|58s|, |; 59s|, — |; 57rf|, |) state. We see that in both 
cases, the distance of the third atom has very little ef- 
fect on the other two atoms: as either Ri or R 2 is 
increased (while keeping the other distance fixed), the 
two stationary atoms consistently demonstrate an attrac- 
tive/repulsive behavior. 

Figure [4] illustrates another type of surface, in which 
there is a significant "ridge" running along one of the 
axes (in this case along the R 2 axis). Such a ridge 
indicates that the two local atoms (e.g. atom 1 and 
atom 2) form a bonded pair, existing even as atom 3 
is moved away. This particular surface corresponds to 




the |59s|,-|;55d|, |;58dl 



asymptotic state. 



FIG. 3. (Color online) Potential energy surface (PES) cor- 
related to the |58s|, |; 59s|, — |; 57d|, |) asymptotic state. 
This surface is analogous to an attractive potential curve for 
the two-body case: As the distance of either the first or last 
atom in the linear chain is increased, the two local atoms re- 
main attracted. The color scheme denotes the energy values 
given in GHz, with the scale presented to the right of the plot. 



B. Potential Wells 

Although the surfaces highlighted in the previous 
section are interesting, ultimately we seek surfaces 
that illustrate potential wells, as these indicate bound 
three-atom systems. Upon separately investigating 
three excited 58s rubidium atoms and three excited 
58p rubidium atoms, we found that in the case of the 
three excited 58p atoms, surface plots corresponding 
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FIG. 4. (Color online) Potential energy surface (PES) corre- 
lated to the |59sk, — 4; 55df , f ; 58d|, — |) asymptotic state. 
We note the "ridge" lying along the R2 axis, which indicates 
that atom 1 and atom 2 are bound (see text). The color 
scheme denotes the energy values given in GHz, with the scale 
presented to the right of the plot. 



to various 58p + 58p + 58p asymptotes illustrated such 
(three-dimensional) potential wells. Specifically, wells 
were determined for the following states: 



11) = 


|58pi 


i;58pi, 


|2> = 


|58pi 


~l;58p 


|3) = 


|58pi 


|;58p|, 


|4) = 


|58pi 


-i;58p 


|5) = 


|58pf 


-§;58p 



2 i 58^2 j — 2! 
2 ' 2 ' 58p^ > 

2 ! 58^2: — 2/ 
2> — 2' 58?2 ' 
2 ' 2 



;58pf : 



As an example, in Figure [5] we show the PES and 
the corresponding two-dimensional projection corre- 
lated to the |1) state, which demonstrates a potential 
well approximately 50-100 MHz deep. Due to the 
large equilibrium separations, e.g. R\ e ~ 22 500 ao 
and i?2e ~ 31 000 ao, we label these bound states 
macrotrimers. 

For the three-atom configuration shown in Figure [lja), 
the (non-zero) oscillation modes can be calculated via: 



(ki + k 2 ) ± ^kf-k^ + k 2 



(7) 



where ki are "effective spring constants" that need to 
be calculated and m is the mass of a single rubidium 
atom. In Figure [ljb) and (c), we show the physical de- 
scription of each eigenmode. Panel (b) corresponds to 
the uj\ eigenvalue and shows that the two outer Rydberg 
atoms vibrate in the same direction, opposite to the in- 
ner atom's direction of motion. Panel (c) corresponds 
to the oj-L eigenvalue and shows that the inner Rydberg 
atom is stationary while the outer two atoms oscillate in 
opposing directions. 

To calculate the ki for Equation |7|, we perform poly- 
nomial fits to two-dimensional cross sections of the po- 
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FIG. 5. (Color online) Potential energy surface (top) 
and two-dimensional projection (bottom) correlated to the 
|58p||;58p|, !;58p|,— |) asymptotic molecular state. The 
main feature of these plots is the three-dimensional well, cen- 
tered at Ri e — 22 500 ao and ifee = 31 000 ao, indicating 
that the three Rydberg atoms are bound together in a lin- 
ear chain (see text). The red (vertical) and blue (horizontal) 
lines imposed on the bottom plot are centered at the well 
minimum and indicate the cross-sections for which quadratic 
fits are performed (see Figure [6|. The color scheme denotes 
the energy values given in GHz, with the scale presented to 
the right of each plot. 



tcntial wells along the i?i and i?2 axes at each well's 
minima (see Fig. [6| . The deepest portions of each well 
can be modeled as a simple harmonic oscillator, and it 
is easily shown that the desired fc-value equals the sec- 
ond derivative of these quadratic fits (with respect to the 
nuclear separation in that particular direction), i.e. 



d 2 V 



(8) 



where Ri e is the equilibrium separation along axis Ri. 
Table [I] summarizes the results of the polynomial fitting 
including the fc e ff-values, the goodness-of-fits (r 2 values) 
of the harmonic oscillator potentials, and the potential 
energy range for which the quadratic assumption is valid. 
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TABLE I. Characteristics of the polynomial fitting procedures 
for each trimer state \i) (see text), assumed to be quadratic. 
A two-dimensional cross-section was taken at the minimum 
of each well for both the Ri and R2 axes, where the well is 
centered (Ri e , R2e)- The k c a -values (in N/m) correspond to 
the numeric fitting along each respective axis. The well depth 
indicates the potential energy range for which the quadratic 
assumption is valid and for which the given r 2 goodness-of-fit 
values are appropriate. 



State 


Axis 


R e (xl0 3 a ) 


k cB (N/m) 


r 2 -value 


Depth (MHz) 


|1) 


Ri 


22.50 


8.37 x 10" 11 


0.9879 


17.50 


R2 


31.00 


4.40 x 10" 12 


0.9926 


16.60 


|2> 


Ri 


23.15 


3.85 x 10" 11 


0.9861 


16.10 


R2 


33.10 


1.07 x 10" 12 


0.9781 


77.50 


|3) 


Ri 


22.80 


2.55 x 10' 11 


0.9908 


37.70 


R2 


34.60 


3.02 x 10" 12 


0.9892 


21.30 


|4> 


Ri 


22.85 


4.63 x 10' 11 


0.9985 


29.30 


R 2 


31.50 


5.54 x 10~ 12 


0.9943 


5.80 


15} 


Ri 


22.00 


6.87 x 10" 11 


0.9982 


16.24 


R2 


31.70 


1.72 x 10~ 12 


0.9992 


11.94 
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FIG. 6. (Color online) Two-dimensional cross sections of the 
three dimensional surface shown in Figure [5] taken along the 
Ri (top) and R2 (bottom) axes at the well's minimum. Each 
cross-section is fitted with a harmonic potential centered at 
the minimum. The results of these fits are summarized in 
Table □ 

Based on these properties, we use Equation ^ and the 
familiar E = fnu(v + 1/2) to find the oscillation energies 
in the deepest portions of the highlighted wells. Tabic |TT| 
lists the vibrational energies of the first few bound states 
for each modal frequency, uj + and w_ . We see that in each 
case, the energies defined by the uj + frequency illustrate 
spacings of about 3-6 MHz, which are separated enough 
to be detected through spectroscopic means. The MHz 



energy values correspond to /xs oscillation periods, which 
are rapid enough to allow for several oscillations during 
the lifetimes of these Rydberg atoms (roughly 500 ^s 
for n ~ 60 [H]). The energies corresponding to the w_ 
frequency are more closely spaced and demonstrate os- 
cillation periods that are slower, but should still be able 
to be detected experimentally. 



TABLE II. Lowest vibrational levels and corresponding bound 
state energies for both oscillation frequencies of each trimer 
state \i) (see text). Energy± represents the bound energies 
associated with ui± and is measured from the bottom of the 
potential well. 



State 


V 


Energy+ (MHz) 


Energy- (MHz) 







2.71 


0.62 




1 


8.13 


1.86 


|1> 


2 


13.55 


3.11 




3 




4.34 




12 
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15.53 







1.84 


0.31 




1 


5.51 


0.92 


|9\ 


2 


9.18 


1.53 


3 


12.86 


2.14 




4 




2.75 




25 




15.60 







1.50 


0.51 




1 


4.48 


1.54 




2 


7.48 


2.57 


|3> 


3 


10.47 


3.60 
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13.47 


4.63 




5 
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16.46 
19.45 


5.66 
6.69 
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0.70 
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2 




2.09 
3.49 
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4.88 







2.46 


0.39 


|5> 


1 
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7.37 


1.17 
1.94 




14 




11.27 
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C. Other asymptotes 



Of course, the potential wells described in section MB 
are not the only wells that exist. During the course of 
our analysis, we also found wells corresponding to vari- 
ous |59s55G?58e?) asymptotes; some examples of which are 
presented below. Although in principle these wells can 
be evaluated in the same amount of detail as was done for 



the |58p58p58p) wells in section MB we merely present 
visual evidence of their existence in this paper. Should 
these additional asymptotes lend themselves to specific 
experimental probing, then computing their respective 
energy levels would be of value. At this time, however, 
such evaluation is beyond the goals of this paper. 
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FIG. 8. (Color online) Potential surface corresponding to the 
1 59s | — |;55<i| — |;58d||) asymptotic state. The color 
scheme denotes the energy values given in GHz, with the scale 
presented to the right of the plot. This particular surface ex- 
hibits two potential wells, each about 30-40 MHz deep. 
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FIG. 7. (Color online) Potential surface corresponding to the 



|59s§ 



|; 55d||; 58d| |) asymptotic state. The color scheme 



denotes the energy values given in GHz, with the scale pre- 
sented to the right of the plot. This particular surface actually 
exhibits a few potential wells, the largest one being between 
40-50 MHz deep. 
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IV. CONCLUSIONS 

The work presented in this letter demonstrates results 
for the = 1/2 symmetry of 58p + 58p + 58p rubidium 
Rydberg atoms. During the course of our analysis, we 
also examined the surfaces corresponding to asymptotic 
states near the Rb 58s + 58s + 58s asymptote, but no 
potential wells were found for this case. The formalism 
could also be applied to the 58c? + 58d + 58<i case, but 
this would correspond to a large increase in computation 
time and thus, is outside the scope of the work shown in 
this paper. Due to an increased number of basis states 
for the triple 58d case, we would expect to find additional 
wells, however. 

In addition, the methodology presented here can be 
applied to Rydberg states of other ^-symmetry values as 
well as to other alkali elements. The current literature 
regarding ultracold multi-body Rydberg physics involves 



FIG. 9. (Color online) Potential surface corresponding to the 
|59s| - §; 55d||; 58d| - \) asymptotic state. The color 
scheme denotes the energy values given in GHz, with the scale 
presented to the right of the plot. The well exhibited in this 
surface is 60-100 MHz deep. 



one Rydberg atom interacting with multiple ground state 
atoms or a ground state molecule. 

We seek to continue to analyze cases of Rydberg 
trimers for various alkali elements, asymptotes, and lin- 
ear Q symmetries, as well as exploring the energy levels 
corresponding to transverse (bending) modes of oscilla- 
tion. Similarly, we seek to extend the theory to differ- 
ent molecular configurations, such as triangular systems. 
The detection of such trimer states could have applica- 
tions in a variety of research areas, including quantum 
information processing and exotic, ultracold chemistry. 

Furthermore, it might be possible to extend the lin- 
ear chain to include 7V-Rydberg atoms, although this 
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is purely speculative at this stage. Such investigations 
could prove fruitful in the advancement of ultracold 
multibody physics. 
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